Identification of rolling resistance as a shape parameter in sheared granular media 
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Using contact dynamics simulations, we compare the effect of rolling resistance at the contacts 
in granular systems composed of disks with the effect of angularity in granular systems composed 
of regular polygonal particles. In simple shear conditions, we consider four aspects of the mechan- 
ical behavior of these systems in the steady state: shear strength, solid fraction, force and fabric 
anisotropies, and probability distribution of contact forces. Our main finding is that, based on the 
energy dissipation associated with relative rotation between two particles in contact, the effect of 
rolling resistance can explicitly be identified with that of the number of sides in a regular polygonal 
particle. This finding supports the use of rolling resistance as a shape parameter accounting for 
particle angularity and shows unambiguously that one of the main influencing factors behind the 
mechanical behavior of granular systems composed of noncircular particles is the partial hindrance 
of rotations as a result of angular particle shape. 
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Most numerical studies on the mechanical behavior of 
granular materials deal with model systems composed of 
disks in 2D or spheres in 3D. This is usually due to the 
technical difficulties that arise when dealing with par- 
ticles of complex shapes in experiments or discrete el- 
ement methods. However, real granular materials are 
rarely composed of spherical particles, and it has been 
shown that the nonspherical shape of the grains strongly 
influences the mechanical behavior of granular systems. 
This influence can be evidenced when characterizing the 
shear strength [H-Q and solid fraction 0-0| : as well as mi- 
crostructural properties such as the distribution of con- 
tact forces [1, Q. The effect of grain shape is thus a 
crucial aspect to be taken into account for a realistic de- 
scription of granular systems. 

One of the numerical "tricks" that can be used to ob- 
tain realistic values of strength and solid fraction while 
using only circular particles in simulations is to partially 
restrict the relative rotations between grains For ex- 
ample, several studies have shown that rolling resistance 
leads to shear strengths and solid fractions that are com- 
parable to those observed in granular soils and rocks. 
However, the extent to which rolling resis- 
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tance can actually be compared to angular shape in more 
general terms, or whether rolling resistance and angular 
shape lead to similar structures at the mesoscopic scale, 
are interesting issues that remain poorly understood. 

In this Letter, we compare, by means of discrete el- 
ement simulations, the effects of rolling resistance and 
angularity. We construct two sets of polydisperse 2D 
packings. In the first set, the packings are composed of 
disks with an increasing magnitude of rolling resistance, 
whereas in the second set, the packings are composed of 
regular polygonal particles of increasing number of sides. 
By comparing various properties extracted from the two 



sets, we find a remarkable matching of the data from 
the disk packings with those of the polygon packings for 
a rolling resistance expressed by a simple equation as a 
function of the number of sides. This one-to-one map- 
ping between the two sets is based on energy dissipation 
considerations and might be generalized to other particle 
shapes. 

All packings are made up of 7500 grains with diameters 
uniformly distributed by volume fractions between 0.6d 
and 2.4d, where d is the mean diameter. In all simula- 
tions, the coefficient of sliding friction fi^ between parti- 
cles is 0.4 and collisions are perfectly inelastic. The par- 
ticles are initially placed in a semiperiodic box lOOd wide 
using a geometrical procedure [l5|. Next, the packing is 
sheared by imposing a constant shear velocity and a con- 
stant confining stress. To avoid strain localization at the 
boundaries, sliding and rolling are inhibited for the par- 
ticles in contact with the walls. The samples are sheared 
up to a large cumulative shear strain 7 = Ax/h = 5, 
where Ax is the horizontal displacement of the upper 
wall and h is the thickness of the sample. All measures 
are averaged over the last 50% of cumulative shear strain 
in order to guarantee that they characterize the behav- 
ior of the system in the steady state, also known as the 
"critical state" in soil mechanics. In all tests, the gravity 
is set to zero. 

The simulations were carried out by means of the con- 
tact dynamics method iB"!^, which assumes perfectly 
rigid particles interacting through mutual exclusion and 
Coulomb friction. For specific implementation of the con- 
tact dynamics method see 

In the first set of samples, composed of disks, the 
rolling_^resistance is introduced through a rolling friction 
law [20[, analogous to the sliding friction law. Although 
rolling friction is introduced here as a numerical param- 
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eter, it may reflect various raaterial parameters such as 
hysteresis, micro-sliding when the elastic moduli are dif- 
ferent, inelasticity (in particular for polymers) and sur- 
face roughness [21j. This law assumes that a contact 
can transmit a torque M not exceeding a limit value 
Mmax = IJ'r^fn, where IS the coefficient of rolling fric- 
tion, £ is the magnitude of the branch vector joining the 
centers of the contacting particles, and /„ is the normal 
force. The scaling of Mmax with i is meant to make fir 
dimensionless. Relative rotation between two grains in 
contact is allowed only if M = Mmax- 

In the second set of samples, composed of regular 
polygonal particles, two types of contact may occur: (1) 
between a corner and a side, and (2) between two sides. 
Side/side interactions represent two constraints and are 
treated by associating two contact points along the com- 
mon side and applying the volume exclusion and the slid- 
ing friction law to each of them. Thus, in practice, two 
contact forces are calculated at each side/side contact. 
However, only their resultant and application point are 
physically relevant, and the result is independent of the 
choice of the two contact points 2^, l23 |. 

The stress components can be calculated from the sim- 
ulation data by the relation cr.y = ndfi^'j), where ric is 
the number of contacts per unit volume and the average 
runs over the contacts c with contact force and branch 
vector [^l- The mean stress ia p — (cri -f cr2)/2, where 
ai and (T2 are the principal stress values, and the devi- 
atoric stress is q = (cti — cr2)/2. It is worth noting that 
in the presence of rolling resistance the stress tensor can 
be asymmetric and a couple-stress tensor may be added 
to the description. However, in all our tests the asym- 
metry is negligibly smaU (i.e., |(cti2 -o'2i)/(o-i2 -l-cr2i)| < 
0.0002), suggesting that we do not need to consider the 
couple stress in the present study. Similar observations 
on the contribution of the couple stress tensor are re- 
ported in ^isj. 

Figure [T] shows the shear strength q/p and solid frac- 
tion v = Vp/V, where Vp is the volume occupied by the 
particles and V is the total volume, as functions of fir for 
the disks and of l/ug for the polygons, where rig is the 
number of sides of the polygons. It can be seen that both 
q/p and v follow similar trends in the two sets as fir and 
l/rig increase. However, a direct comparison of the data 
between the two sets is not possible in this representation 
due to the different physical meanings of fir and l/n^. 

The respective effects of rolling friction and angular 
shape can be compared by their roles in the hindering of 
relative rotation. Let us consider a particle (a disk with 
rolling friction and a regular pentagon) that rolls on a 
horizontal plane with a vertical force N exerted at its 
center of mass and that is not allowed to slide; see Fig. 
[SJa). Figure [2l[b) shows the horizontal force T that must 
be applied at the center of mass in order to make the 
particle roll, as a function of the rotation angle 6. The 
work needed to displace the particle a distance equal to 
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FIG. 1: (color online) Shear strength q/p (Up) and solid frac- 
tion f (Down) as functions of fir for the disks and of 1 /us for 
the polygons. Error bars indicate the standard deviation. 




FIG. 2: (a) Schema of rolling on a plane, (b) Horizontal force 
T required for rolling a distance equal to the perimeter, (c) 
Trajectory of the center of mass of the polygon (dashed line) 
and definition of the mean dilatancy angle tp. 



its perimeter is 



Wd = AlTflrRdN 



(1) 



for the disk with rolling friction, where Rd is the radius 
of the disk and the magnitude of the branch vector £ 
(necessary to calculate Mmax) has been replaced by the 
disk diameter, and 



Wp — ns{l — cos{n /ns))RpN 



(2) 



for the polygon, where Rp is the radius of its circumcircle. 
Assuming equal work, i.e. Wd = Wp, we arrive at the 
following mapping between fir and Us 



fir — (1/4) tan V', 



(3) 



where it has been assumed that both particles have the 
same perimeter (i.e., Rp = Rdijr/ns)/ siii{'iT/ns)), and 
ijj = 7r/(2ns) is the mean dilatancy angle of the trajectory 
of the center of mass of the polygon (see Fig. [UJc)). For 
a similar attempt to quantify the role of grain shapes in 
hindering relative rotation, see (26| . 

Figure 13] shows the shear strength q/p and solid frac- 
tion v as functions of fir for the disks and of (1/4) tan?/' 
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FIG. 3: (color online) (Up) Shear strength q/p as a function 
of fir for the disks and of (l/4)tan'(/) for the polygons, both 
from raw simulation data (full symbols) and as predicted by 
Eq. 13] (empty symbols). (Down) Solid fraction i/ as a function 
of fj.r for the disks and of (1/4) tan-i/i for the polygons. Error 
bars indicate the standard deviation. 




for the polygons. Remarkably, the shear strengths and 
solid fractions of the two sets of packings collapse, both 
increasing and decreasing, respectively, with /i^ and 
(1/4) tan?/; and tending to a constant value at fir — 
(1/4) tan?/' ~ 0.1. In other words, from a macro-scale 
viewpoint, a packing of regular polygons of Ug sides is 
equivalent to a packing of disks with a coefficient of 
rolling friction /i,. given by Eq. [3) This result supports 
also the choice of the required energy for rolling as a 
relevant physical quantity for the rheology of granular 
materials. 

The mapping evidenced in Fig. [3] hints at similar pack- 
ing structures in the two sets. Figure |3] shows two snap- 
shots: one representing a disk packing with fir — 0.05 
and the other representing a polygon packing with Hs ~ 8 
(note that 0.05 ~ (1/4) tan(7r/(2 * 8)) according to Eq. 
[3|). The contact forces are represented by segments join- 
ing the particle centers, with a thickness proportional to 
the force magnitude. We observe that the force-carrying 
backbone is astonishingly similar in the two systems. 

From the expression of the stress tensor, it can be 
shown that the shear stress q/p reflects the packin g st ruc- 
ture and force transmission via a simple relation j27| : 



q/p ~ (l/2)(ac + a„ + at), 



(4) 



where Oc, a„, and at, are the anisotropics of the angular 
distributions of contact orientations Pn{d), normal forces 
{fn){d), and tangential forces {ft) (9), respectively, as a 
function of contact orientation 9, which are approximated 
by their lowest order Fourier expansions: 



Pni9) ~ l/Tr{l + a,cos2{9~e,)}, 
{fn){9) ^ (/„){l + a„cos2(0-0„)}, 
{ft){9) ^ -{fn)atsm{0~Ot), 



(5) 



where (/„) is the mean normal fore, and 9c = 9n = 9t 
are the corresponding privileged directions, which, in the 
steady state, coincide with the principal stress direction. 



FIG. 4: Snapshots of the force network in (a) a system com- 
posed of disks with rolling friction {fir = 0.05) and (b) a 
system composed of octagonal particles. The line thickness is 
proportional to the normal force. The floating particles are 
represented in light grey. 
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FIG. 5: (color online) Contact anisotropy ac (Up) and 
force anisotropics, a,i (full symbols) and at (empty symbols) 
(Down), as functions of fir for the disks and of (1/4) tamp for 
the polygons. Error bars indicate the standard deviation. 



Equation 3] reveals distinct origins of the shear strength 
in terms of force and texture anisotropy. The empty sym- 
bols in Fig. [3] show q/p as predicted by Eq. |31 We see 
that this equation approximates well the shear strength 
for all raw data. 

The anisotropics a^, an, and at are shown in Fig. [5] 
as functions of fir for the disks and of (1/4) tan -0 for 
the polygons. It is remarkable that all anisotropics are 
almost identical between the two sets. This correspon- 
dence is only broken for polygons with small numbers 
of sides, i.e., for Ug = 3 and 4. This happens because 
for these polygons the contact orientation is strongly in- 
fluenced by the low rotational symmetry of the particles 
and the orientations of the sides rather than the relative 
positions of the particles. 

The mapping between rolling friction and angular 
shape of particles is also reflected by the probability den- 
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FIG. 6: (color online) Probability distribution function of nor- 
malized normal forces fn/{fn} for the two systems of Fig. [4] 
The inset shows the evolution of the exponents a and P (see 
text for definitions) with /ir for the disks and with (1/4) tan-i/) 
for the polygons. 

sity function (PDF) of normal forces displayed in Fig. 
[6l In this figure, we compare the PDFs of the two sam- 
ples shown in Fig. The two PDFs are almost identi- 
cal. The PDF can be approximated by a power law (i.e., 
PDF cx (/„/ (/n))~") in the range of small forces and by 
an exponential function ( i.e., PDF oc e^^^~^"/^^"'^^) in 
the range of strong forces [28| . The inset in Fig. [6] shows 
the coefficients a and /? as functions of jir for the disks 
and of (l/4)tan?/i for the polygons, confirming that the 
similarity of the PDFs is maintained for the whole range 
of rolling frictions and numbers of sides studied in this 
work. 

To sum up, the simulations presented in this Letter 
provide strong evidence for the mapping between the two 
studied parameters, i.e., rolling resistance and shape an- 
gularity. This correspondence was established by con- 
sidering shear strength, solid fraction, force and fabric 
anisotropics, and the PDFs of normal forces in the steady 
state. A practical consequence of this finding is that 
rolling resistance may be employed to imitate the effect 
of angular shape in discrete-element simulations of gran- 
ular materials. More importantly, it suggests that the 
hindrance of particle rotations is a major effect of angu- 
lar particle shape. In this picture, the effect of rolling 
friction is to force the particles to rearrange as if they 
were glued to each other. In this way, the clusters of two 
or more particles behave as non-circular particles. This 
result may be tested in other loading conditions and it is 
potentially extensible to other particle shapes (in 2D and 
3D), opening new scopes in modeling complex granular 
systems. 
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